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ABSTRACT 


In  this  third  part  of  A  Theory  of  Spectral  Line  Shape  for  Atmospheric 
Absorption,  an  evaluation  of  the  frequency  dependence  of  the  complete  Fano 
collision  operator  has  been  made  using  the  collision  time  theory  developed 
earlier.  This  results  in  a  considerable  reduction  of  the  terms  in  the  Fano 
operator.  An  analysis  of  the  spectral  line  wing  shows  that  there  are  two 
different  regimes:  the  near  wing  and  the  far  wing.  These  regimes  are 
demarcated  by  a  turning  point  at  a  frequency  displacement  of  kT/h.  In  the 
very  far  wing,  beyond  the  turning  point,  the  behaviour  of  the  Fano  operator 
changes  and  the  effects  of  binary  collision  overlap  and  photon  absorption 
into  the  binary  collision  complex  become  important.  The  line  wing  behaviour 
is  corroborated  qualitatively  by  experimental  results  for  carbon  dioxide  near 
4 . 2  /Lim. 


r£sum£ 


Dans  cette  troisieme  partie  de  notre  theorie  des  profils  de  raies 
spectrales,  une  evaluation  de  la  dependance  par  rapport  a  la  frequence,  de 
l'operateur  de  collision  de  Fano,  dans  sa  totalite,  est  elaboree  grace  a 
1 ' application  de  la  theorie  de  la  duree  de  collision  d^veloppee 
ant^rieurement .  Le  r6sultat  donne  une  importante  reduction  des  termes  qui 
composent  l'operateur  de  Fano.  Une  analyse  de  l'aile  de  la  raie  spectrale 
demontre  qu'il  y  a  deux  regimes  differents:  l'aile  proche  du  centre  et  l'aile 
lointaine.  La  demarcation  entre  les  deux  regimes  s'identifie  par  un  point 
tournant  qui  se  situe  a  un  deplacement  en  frequence  egal  k  kT/h.  Dans  l'aile 
lointaine,  au-delA  du  point  tournant,  le  comportement  de  l'operateur  de  Fano 
change  et  des  effets  de  chevauchement  de  collisions  binaires  ainsi  que 
1' absorption  d'un  photon  par  le  complexe  de  la  collision  binaire  deviennent 
importants .  Le  comportement  de  l'aile  de  la  raie  est  confirme  de  fagon 
qualitative  par  des  resultats  experimentaux  pour  le  bioxide  de  carbone  dans 
la  region  des  4.2  Jim. 
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EXECUTIVE  SUMMARY 


The  performance  of  electro-optical  sensing  devices  in  both  the  military 
and  the  civilian  remote  sensing  fields  can  be  limited  by  the  propagation 
characteristics  of  the  atmosphere.  It  is  often  assumed  that  the  attenuation 
by  molecular  absorption  is  sufficiently  well  known  for  it  to  be  completely 
predictable.  Unfortunately,  this  is  not  the  case.  The  two  main  atmospheric 
transmission  models,  LOWTRAN  (for  broad-band  applications)  and  FASC0D2  (for 
narrow  band  and  lasers),  still  have  certain  limitations. 

Recent  transmission  studies  at  DREV  have  revealed  that  both  models 
deviate  from  measured  transmission  in  the  two  main  atmospheric  transmission 
windows.  One  feature  which  is  still  modelled  inadequately  is  the  absolute 
transmission  of  the  water  vapour  continuum  over  the  wide  range  of 
temperatures  in  the  troposphere.  Another  one  is  the  combined  carbon  dioxide 
and  nitrogen  continuum  in  the  4.2  juun  region.  A  third  problem  concerns  the 
edges  of  strong  bands  of  water  vapour. 

Two  of  the  most  important  absorption  line  phenomena,  which  have  not  yet 
been  completely  described  theoretically,  are  line  coupling,  or  overlap,  and 
the  shape  of  line  wings  far  from  the  line  centres.  These  phenomena  are 
especially  relevant  to  the  above  problem  areas  of  LOWTRAN  and  FASC0D2 . 

In  this  third  part  of  the  theory  of  spectral  line  shapes  for  infrared 
atmospheric  transmission  the  problem  of  the  complete  wing  of  an  isolated  line 
is  examined.  The  collision  time  theory,  developed  in  Part  II  to  describe  the 
frequency  dependence  of  line  widths  and  shifts,  is  applied  to  the  complete 
Fano  operator.  Two  distinct  regions,  the  near  and  far  wings,  having  different 
frequency  behaviour  are  identified.  The  theoretical  line  shape  is  in  accord 
with  experimental  results  for  the  continuum  absorption  of  carbon  dioxide  near 
4.3  jum. 

This  work  is  part  of  a  programme  to  improve  transmission  prediction, 
involving  experimental  and  theoretical  studies  at  DREV  and  modelling  at  the 
Centre  for  Research  in  Earth  and  Space  Science,  York  University. 
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1.0  INTRODUCTION 

The  non-Lorentzian  nature  of  the  far  wings  of  spectral  line  shapes  is  an 

important  phenomenon  in  the  description  and  modelling  of  atmospheric 

transmission.  The  CO  continuum  near  4.3  /xm  and  the  H  O  vapour  continuum, 

2  2 

which  pervades  the  atmospheric  transmission  windows  in  the  3-5  and  8-12  /xm 
regions,  are  good  examples  of  the  cumulative  effects  of  non-Lorentzian  line 
wings.  Fano  (Ref.  1)  developed  an  approach  to  spectral  line  shape  theory 
whose  essential  basis  was  the  definition  of  a  collision  super  operator  m{0)) 
which  is  a  function  of  the  frequency  of  the  applied  radiation  field.  The 
diagonal  and  non-diagonal  matrix  elements  of  m(ct))  produce  frequency- dependent 
line  widths  (and  shifts)  and  coupling  parameters,  respectively.  Various 
approaches  have  been  used  to  describe  the  non-Lorentzian  nature  of  the  far 
wings  of  lines  (Refs.  2-10).  Few  of  these  treat  directly  the  composition  of 
the  m(o))  operator  as  a  function  of  frequency.  In  some  cases,  empirical  models 
are  used  to  introduce  a  collision  time  as  a  parameter  in  the  far  wing 
behavior.  An  exception  is  Rosenkranz ' s  theory  (Ref.  8)  and  extensions  of  his 
theory  (Refs.  9,10).  This  is  a  statistical  theory  in  which  it  is  assumed  that 
the  molecules  behave  as  though  they  were  stationary  in  the  region  far  from 
the  line  centre:  the  collisions  are  assumed  to  have  infinite  duration.  In 
addition,  the  distance  from  the  line  centre  is  taken  as  much  greater  than  the 
distance  between  individual  lines.  Because  of  these  assumptions  this  theory 
applies  asymptotically  to  the  extreme  far  line  wing.  Furthermore,  the 
approximation  for  the  Fano  operator  results  in  a  band-averaged  relaxation 
parameter  which  takes  into  account  line  coupling  without  the  necessity  of 
calculating  the  complete  set  of  coupling  coefficients. 

Until  recently,  little  if  any  work  has  been  done  on  the  detailed 
analysis  of  the  frequency  dependence  of  the  complete  Fano  operator.  Neither 
has  any  attempt  been  made  to  perform  line  shape  calculations  without  first 
making  quite  drastic  approximations  to  m(u))  .  This  is  because,  apart  from  the 
fact  that  m(w)  is  composed  of  a  formidable  number  of  terms,  there  existed  no 
simple  methods  for  explicating  the  frequency  dependence  of  the  operator, 
which  is  manifested  primarily  in  terms  of  of f- the -energy  shell  transition 
operators . 
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In  two  recent  papers  (Refs.  11,12)  we  have  developed  an  approach  to 
spectral  line  shape  applicable  to  the  conditions  for  atmospheric  infrared 
transmission.  The  theory  was  formulated  in  Ref.  11  (henceforth  referred  to  as 
Part  I)  and  was  based  on  a  kinetic  equation  for  the  density  matrix  for  a 
system  of  absorbing  molecules  in  a  thermal  bath.  An  iterative  solution  for 
the  linearized  kinetic  equation  was  developed  as  a  density  expansion  about  an 
isolated  line  approximation  based  on  Fano's  collision  operator.  The 
absorption  coefficient  was  determined  using  linear  response  theory,  together 
with  the  fluctuation-dissipation  theorem  and  detailed  balance.  The  isolated 
line  shape  has  non-Lorentzian  wings  due  to  frequency- dependent  shift  and 
width  parameters,  given  by  the  real  and  imaginary  parts,  respectively,  of  the 
diagonal  elements  of  an  averaged  m{u))  .  New  expressions  were  found  for  the 
line  coupling  terms,  dependent  on  the  line  strengths  of  coupled  lines  as  well 
as  the  non-diagonal  element  of  m(u))  . 

In  the  second  paper  (Ref.  12)  (henceforth  referred  to  as  Part  II), 
formal  scattering  theory,  in  the  frequency  domain,  was  used  to  derive  a 
general  relation  between  the  of f- the -energy  shell  transition  operator  (which 
forms  the  primary  basis  of  m(to)  )  and  its  on-the-energy  shell  counterpart. 
This  relation .  involves  a  spectral  distribution  function  $(E,A6t))  composed  of  a 
complex  collision  time  operator  T(E)  ,  a  time  operator  2T  (E)  associated  with 
the  free  propagator  (E-H  )  ,  and  Aw  the  displacement  from  the  line  centre. 

O 

In  Part  II,  this  theory  was  applied  only  to  the  near  line  wing  where 
displacements  from  the  line  centre  are  small  compared  with  the  average 
kinetic  energy  kT.  This  condition  allowed  the  retention  of  only  those  terms 
in  m(o))  which  are  related  to  the  usual  impact  approximation  parameters  such 
as  the  shift  A(0)  and  width  T(0). 

The  main  objective  of  this  work  is  to  use  the  collision  time  theory 
developed  in  Part  II  to  carry  out  an  analysis  of  the  contributions  of  all  the 
terms  in  the  Fano  operator  to  the  spectral  line  wing,  including  the  very  far 
wing  region  not  considered  in  Part  II.  We  start  with  a  review  of  the  nature 
of  the  terms  in  the  Fano  operator  as  discussed  by  that  author.  This  is 
followed  by  a  summary  of  the  collision  time  theory  results  of  Part  II.  We 
then  apply  this  theory  to  each  of  the  terms  of  the  Fano  operator  and  examine 
the  contributions  to  the  near  and  far  line  wing. 
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This  work  was  carried  out  at  DREV  between  November  1993  and  June  1994 
under  PSC  32E.  EO  Battlefield  Surveillance. 

2 . 0  THE  COMPLETE  FANO  OPERATOR 

The  original  Fano  operator  m(w)  was  given  by  eq.  49  of  Ref.  1. 

m(w)  =  t(w+H*)  -  t*  (H  -w) 

o  o 


27ri 


r»+i  V 
dip 
oo+  i->7 


l/f-H  ip-0)-E 


t{ip)  t  {ip-u) 


ip-H.  ip-u>-E 


,  [1] 


where  H  is  the  interaction-free  Hamiltonian  and  t(E)  the  transition 
o 

operator.  In  Parts  I  and  II  we  took  translational  motion  into  account  and  in 

this  case  we  have  a  thermally  averaged  operator  M(Aw)  ,  where  the  average  is 

taken  over  the  thermal  bath  molecules  which  are  assumed  to  have  no 

interaction  with  the  applied  radiation.  In  the  frequency  domain  the  variables 

are  Ao)  ~  (u>  -  a>  )  :  the  frequency  displacement  from  a  line  centred  at  a> 

ab  ab  ab 

After  integration  over  \j)  (Ref.  1),  the  matrix  elements  of  M(Aw) ,  are 
represented  by 


where 


M ( Aw)  =  <  m( Au>)  >  —  <  F  + 

acdb  acdb  AV 


I  dk'  {S  +  C  +  P  +  C  +  P  } 

J  1  1  2  2J 


AV 


F  =  t  (e  +E  +hAw  )  -  t  (e  +E  -hAu>  )  , 

aa  a  k  ab  bb  b  k  ab 


[2] 

[3] 


S  =  +iir  6  [E  -E  -hAw  ]  t  (e  +E  )  t  (e  +E  ) 
k#  k  cb  ac  c  k'  db  b  k 


+i?r  S[ E  -E  +hAw  ]  t  (e  +E  )  t  (e  +E  ) 

k'  k  ad  ac  a  k  db  d  k' 


[4] 


C1  =  +  2 


t  (e  +E  )  t  (e  +E  -hAw  )-t  (e  +E  )  fc  (e  +E  -hAa>  ) 
ac  c  k'  db  d  k'  cd  ac  a  k  db  b  k  ab 

E  -E  -  h0) 
k '  k  ac 


t  (e  +E  +hA 0)  ) t  (e  +E  )-t  (e  +  E  +hAto  )t  (c  +E  )  1 
ac  c  k'  cd  db  d  k '  ac  a  k  ab  db  b  k 

E  -E  +  hCt) 
k '  k  db 


[5] 
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In  the  preceding  paper.  Part  II,  a  reduced  Fano  operator  was  used 
consisting  only  of 

M(Aw)  “  <  F  +  I  dk'S  > 

acdb  J  AV 

This  form  is  valid  in  the  near  wing  when  Aw  is  small.  As  pointed  out  by  Fano 

(Ref.  1),  the  remaining  terms  C^,  C^,  and  P ^  depend  on  transition  matrix 

elements  evaluated  of f- the -energy  shell.  In  the  above  form,  the  off-shell 

excursions  are  composed  of  energy  increments  ±hAw  and  ±hAct>  .  When  these 

ab  cd 

are  negligibly  small,  as  in  the  impact  approximation  and  in  the  near  line 
wing,  these  terms  vanish.  Fano  attributed  these  terms  to  the  effects  of 
transient  stages  of  collisions,  as  opposed  to  the  terms  F  and  S  which  relate 
and  are  twice  Fano 7 s  due  to  a  misprint  in  Fano's  eq.  55a. 
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to  completed  or  well-separated  collisions.  The  implication  is  that  the 
transient  effects  arise  from  long  collision  times  so  that  a  collision  is  not 
completed  before  the  start  of  another  one.  We  can  expect  these  transient 
effects  terms  to  become  important  in  the  very  far  line  wing.  A  major  task  in 
this  paper  is  the  determination  of  these  terms  using  the  collision  time 
theory  developed  in  Part  II. 

3.0  SUMMARY  OF  COLLISION  TIME  THEORY 


In  Part  II  we  derived,  through  formal  scattering  theory,  operator 
expressions  relating  of f- the- energy  shell  transition  operators  to  their 
on-shell  values.  Thus  we  showed  that 

t (E+AE+iO )  =  t(E+iO)  $(E,AE),  [9] 

where  the  spectral  distribution  operator  $(E,AE)  is  given  in  terms  of  two 
time  operators  T(E)  and  J (E)  by  the  relation 


$(E, AE) 


=  |l  -  iAE  J(E)/h  j  |l  - 


iAE  [T(E)  +9'<E)  ]  /h 


f- 


[10] 


The  argument  AE  corresponds  to  the  energy  increment  of f- the- energy  shell.  In 
the  spectral  line  application  AE  =  hAa),  and  the  on-shell  energy  will 
correspond  to  half-shell  arguments  as  in  eq.  2.  The  complex  collision  time 
operator  T(E)  is  defined  by 


ih 


T[E)  =  Lt 


t  (E+iT)) 


[11] 


T)-*0  (E-  H 


and  the  other  time  operator  J  (E)  is  associated  with  the  free  pair  propagator 


(E-H  )  by  the  expression 
o 


ih 


7 (E)  =  Lt 


[12] 


(E-  H  +i7)) 
o 


Another  convenient  form  for  $  is 

$(E,AE)  =  1  +  (iAE/h)  T(E) 


1  -  i(AE/h)T(E) 


[10a] 
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where 


T(E)  =  T(E)  +  J  (E) 


The  function  $(E,AE)  is  a  spectral  distribution  function  for  the  energy 
increment  AE  which  compares  the  total  time  T(E)  the  molecular  pair  spend  in 
each  other's  vicinity  during  collision  with  the  time  9"  (E)  they  would  have 
spent  in  the  interaction  vicinity  in  the  absence  of  the  collisional 
interaction.  The  asymptotic  limit  of  $(E,AE),  as  AE  tends  to  infinity,  is 


Lt  $(E,AE)  =  J  (E)  [T(E}+J(E)] 

Ae-*oo 


[10b] 


4.0  APPLICATION  OF  COLLISION  TIME  THEORY  TO  THE  COMPLETE  FANO  OPERATOR 


In  order  to  evaluate  M(Aw)  we  will  make  one  simplifying  assumption.  In 
accordance  with  the  first  order  perturbation  solution  for  $  of  Part  II,  only 
the  diagonal  matrix  elements  of  the  operators  T(E)  and  JT{E)  are  retained. 
Also,  since  they  will  always  be  on  an  energy  shell,  we  will  drop  the  E 
argument  and  merely  write  $(AE)  . 


We  are  now  in  a  position  to  calculate  all  the  terms  of  eq.  2.  Using  eq. 
9  the  forward  scattering  terms  F  (eq.  3)  are  given  by 

F  =  t  (e  +E  +hAw  )  -  t  (e  +E  -hAw  ) 
aa  a  k  ab  bb  b  k  ab 

=  t  (e  +E  )$  (Aw  )  -  fc  (e  +E  )  $  (-A w  )  [14] 

aa  a  k  a  ab  bb  b  k  b  ab 

Similarly,  the  singular  terms  S  of  eq.  4  can  be  written 

S  =  +i n  8  [E  -E  -hAw  ]  t  (e  +E  )  t  J (e  +E  ) 
k'  k  cb  ac  c  k'  db  b  k 

+i7t  5  [E  -E  +hAw  ]  t  (e  +E  )  t*(e  +E  ) 

k'  k  ad  ac  a  k  db  d  k' 

=  +iirr  8  [E  -E  -hAw  ]  t  (e  +E  +hAw  )  t  (e  +E  ) 

k'  k  cb  ac  a  k  ab  db  b  k 

+i7T  8[ E  -E  +hAw  J  t  (e  +E  )  t*ie+E-hAu>) 

k'  k  ad  ac  a  k  db  b  k  ab 


=  in  8[E  -E  -hAw  1®  (Aw  J 

k '  k  c  b  a  ab 


+  8[ E  -E  +hAw  ]$  (“Aw  )  t  (e  +E  )  fc  (e+E  )  . 

k'  k  ad  b  ablacakdbbk 
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Next,  the  C  terms  of  eq.  5  can  be  expressed  as  follows: 


ci  - 1 


t  (s  +E  )  t  (e  +E  +hw  -hw)-t  (e  +E  )  t  (e  +E  -hAw) 

ac  c  k'  db  d  k'  cd  ac  a  k  db  b  k 

E  -E  -hw 
k '  k  ac 


t  (e  +E  +hw-hw  )  t  (e  +E  )-t  (e  +E  +hAw) t  (e  +E  )  1 

ac  c  k#  cd  db  d  k 7  ac  a  k  db  b  k 

E  -E  +hW 
k'  k  db 


1 1 

2 1  '.ivv'yvv 

$*(-Aw  ) 

b  ab 

$  (Aw  ) 

a  ab 

E  -E  -hw 
-  k'  k  ac 

E  -E  +hw 
k '  k  db  J 

-  t  (e  +E  ) t  (c  +E  ) 

ac  c  k'  db  d  k' 

$  (-Aw  ) 

d  cd 

$  (Aw  ) 

c  c  d 

E  -E  -hw 
*-  k '  k  ac 

E  -E  +hW 
k '  k  db  J 

first  set  of  principal 

value  terms 

P  (eq.  6)  v 

[16] 


off-shell  t  elements  in  terms  of  on-shell  elements  (using  eq.  9) : 


P  =  -  i  P 
1  2 


t  (c  +E  )  tJe  +E  )®  (-Ao>  )  +  t  (€  +E  )  t  (e  +E  )®  (Aw  ) 

acc  k'  db  d  k'  d  cd  aca  k  db  b  k  a  ab 

E  -E  -hAw 
k'  k  cb 


-  P 


t  (e+E  )t  (e  +E  )$  (Aw  )+t  (e  +E  )  t  (e  +E  )$  (-Aw  ) 

ac  c  k'  db  d  k'  c  cd  ac  a  k  db  b  k  b  ab 

E  -E  +hAw 

k '  k  ad 


_  ,  t  (e  +E  ) t  (e  +E  ) 
21  ac  c  k'  db  d  k' 


$  (-Aw  ) 

d  cd 


E  -E  -hAw 

k f  k  cb 


$  (Aw  ) 

c  cd 

E  -E  +hAw 

k '  k  ad 


t  (e  +E  ) fc  (C  +E  ) 

ac  a  k  db  b  k 


$  (-Aw  ) 

b  ab 


E  -E  +hAw 
k '  k  ad 


$  (Aw  ) 

a  ab 

E  -E  -hAw 
k '  k  cb 


[17] 


In  the  last  set  of  terms,  (eq.  7)  and  P ^  (eq.  8),  we  first  use  eq.  9 
to  write  the  t  matrices  in  terms  of  their  appropriate  on-shell  counterparts. 
At  this  point  we  can  identify  four  integrals  over  the  variable  \p  and  write 
and  in  the  following  form: 
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C2  =1 


t  (e  +  E  )t(e  +  E  ) 
ac  c  k '  db  d  k ' 


.  E  -E  -hw  E  -E  +hw 

L  k'  k  ac  k'  k  db 


-t(e  +  E)t(e  +  E) 

ac  a  k  db  b  k 


E  -E  -hCt)  E  -E  +hw  , 
•  k'  k  ac  k'  k  dbJ 


P2  2 


fc  (c  +  E  )  t  (G  +  E  ) 
ac  c  k'  db  d  k' 


,  E  -E  -hAw  E  -E  +hAw  . 

L  k'  k  cb  k'  k  adJ 


-  fc  (G  +  E  )  t  (G  +  E  ) 
ac  a  k  db  b  k 


E  -E  -hAw  E  -E  +hAw 
L  k'  k  cb  k'  k  ad* 


where  the  four  integrals  1^,  I^,  1^  and  1^  are  all  of  the  same  type 

evaluated  in  Appendix  A: 

oo 

r- 


i  =4 

1  7T1 


dip  - 

0-G  -E 

-00  ^  c  k' 


$  -E  )  $  (lp-£  -E,  -hAw  ) 

c  c  k'  d  d  k'  cd 


=  2©  ^Au>  J  +  <M-AC(>  J  , 


cd  cd  d  cd 

„00 


1=4 

2  7T1 


d \p 


.  \p-c  -E 

^00  r  a  k 


$  (1 Jj-£  -E  )  $  (0-£  -E  -hAw  ) 

a  a  k  b  a  k  ah 


I  = 


_1 

TTi 


=  20  (Aw  )  +  $  (-Aw  )  , 

ca  ab  b  ab 

P 


,oo 


di j)  - 

t b-e  -E  -hAw 

00  c  k '  cd 


<E>  ( il>-e  -E  )  $  (i/j-e-E  -hAw  ) 

c  c  k'  d  d  k'  cd 


[2$  (Aw  )  +  $  (Aw  )  ]  , 

cd  cd  c  cd 


00 


i=4 

4  ni 


difi 


„  <A-e  -E,  -hAw 

-oo  a  k  ab 


$  (i b-e  -E  )  $  (lA-e  -E  -hAw  ) 

c  a  k  d  ak  ab 


[2$  (AW  )  +  $  (AW  )] 

cd  ab  a  ab 


[18] 

,  [19] 

and  are 

[20a] 

[20b] 

[20c] 

[20d] 


The  new  spectral  distributions  appearing  in  eg.  20  are: 
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B^Aw 


cd 


$  j^Aw 


cd 


(1 +  T  / T  )  (1-T  /T  )/2 

c  c  d  d 


(1-T  /T  )  (1+T*/T* )  /2 
c  c  d  d 


(T  T  /T  ) 
c  d  d 


(  T  +T  )  -i T  T  Aw 
d  c  c  d  cd 


+ 


(T  T*/T  ) 
c  d  c 


—  *  _  _  _  * 

( T  +T  ) -iT  T  Aw 

d  c  c  d  cd 


®*(-Aw  J  [21] 

d  cd 


$  (Aw  ) .  [22] 

c  cd 


It  is  important  to  note  the  different  nature  of  these  functions  compared 

with  $  (Aw  )  and  $  (-Aw  )  .  Both  0  (Aw  )  and  ^  (Aw  )  are  mixed  functions 
c  cd  d  cd  cd  cd  cd  cd 

—  ★  —  ★ 

of  the  collision  times  T  ,  T  ,  T  and  T  :  it  is  not  possible  to  reduce  these 

c  c  d  d 

into  separate  parts  dependent  on  the  collision  time  pairs  T  ,  T  and  T  ,  T 

*  c  c  d  d 

alone,  as  in  $  (Aw  )  and  $  (-Aw  ).  Terms  containing  functions  ©  (Aw  )  and 
c  cd  d  cd  cd  cd 

^  (Aw  )  therefore  refer  to  overlapping  binary  collisions  since  it  is  not 
c  a  ca 

possible  to  distinguish  between  each  binary  collision  through  its  time  of 
duration.  This  is  in  accord  with  Fano's  reference  to  transient  effects  (Ref. 
1)  arising  from  collisions  which  cannot  be  regarded  as  entirely  isolated 
events . 


The  asymptotic  limits  for  the  different  spectral  distribution  functions, 

as  Aw  tends  to  infinity,  are  important: 
cd 

Lt  $  (Aw)  =  (ST  /T  )  ,  [10c] 

•  c  c  c 

Aw->oo 

Lt  $*(-Aw)  =  (J*/T*),  [lOd] 

.  d  d  d 

Aw-»oo 

Lt  0  (Aw)  =  -  (1-T  /T  )  (1-T*/T*)  /2  =  Lt  *  (Aw)  [21a] 

A  cd  c  c  d  d  .  cd 

Aw**  oo  Aw~*  oo 

Substituting  for  eq.  20  in  eqs .  18  and  19,  and  collecting  and  regrouping 

the  terms  C  and  C  (eqs.  16  and  18)  and  P  and  P  (eqs.  17  and  19)  we  can 
12  12 

write 


c  =  c  +  c 
1  2 


1 

2 


-  t  (e  +E  ) t  (e  +E  ) 

ac  c  k '  db  d  k' 


"  $*(-Aw  )-I 

d  cd  1 

,  E  -E  -hw 

L  k#  k  ac 


$  (Aw  J+I 

c  cd  3 

E  -E  +hw 
k'  k  db 


+  t  (e  +E  )  t  (e  +E,  ) 

ac  a  k  db  b  k 


$  (-Aw  )-I 

b  ab  2 

E  -E  -hW 
k '  k  ac 


$  (Aw  ) +1 

a  ab  4 

E  -E  +hw 
k'  k  db 
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t  (e  +E  ) fc  (e  +E  ) 

ac  c  k'  db  d  k' 


0  (Aw  ) 

cd  cd 


#  (Aw  ) 

cd  cd 


E  -E  -hw  E  -E  +hw 
k '  k  ac  k '  k  db 


t  (e  +E  ) t  (e  +E  ) 

ac  a  k  db  b  k 


0  (Aw  )  (Aw  ) 

cd  ab  cd  ab 


E  —E  — hw  E  -E  +hw 

k '  k  ac  k '  k  db 


P  =  P  +  P  =  -  t  (e  +E)  fc  (e  +E1  ) 

1  2  2  I  a  c  c  k '  db  d  k ' 


-#  (Aw  )-l  $  (-Aw  )-I 

c  cd  3  d  cd  1 

-  +  - 

E  -E  +hAw  E  -E  -hAw 

k'  k  ad  k'  k  cb 


fc  (£  +E  ) fc  (e  +E  ) 

ac  a  k  db  b  k 


$  (-Aw  ) -I  $  (Aw  )+I 

b  ab  2  a  ab  4 


E  -E  +hAw 

k '  k  ad 


E  -E  -hAw 

k '  k  cb 


fc  (e  +E  ) fc  (e  +E  ) 

ac  c  k'  db  d  k' 


$  (Aw  ) 

cd  cd 


0  (Aw  ) 

cd  cd 


E  -E  +hAw  E  -E  -hAw 

k'  k  ad  k'  k  cb 


+  fc  (0  +E  ) fc  (e  +E  ) 

ac  a  k  db  b  k 


0  (Aw  ) 

cd  ab 


V  (Aw  ) 

cd  ab 


E  -E  +hAw  E  -E  -hAw 

k '  k  ad  k '  k  cb 


Substituting  in  eq.  2  for  the  terms  F,  S,  C  and  P  from  eqs .  3,  4,  23  and 
24,  respectively,  taking  the  thermal  average  over  the  reservoir  molecules  and 
re-introducing  the  momentum  labelling  for  the  transition  matrices,  we  obtain 
the  matrix  elements  of  M(Aw) 


M  { Aw)  =  <  F  +  |  dk '  {S  +  C  +  P  +  C  +  P}> 
acdb  J  112  2  AV 


=  Tr  dk  feq(k) 
2 

2 


X  <  <k  t  (0  +E  )  k>  $  (Aw  )  6  -<k  t  (e  +E  )  k>  *  (-Aw  )6 

I  aa  a  k  a  ab  db  bb  b  k  b  ab  ca 


+  dk'  ni  <k  fc  (e  +E  )  k'xk'  t  (e  +E  )  k> 

ac  a  k  1  db  b  k  1 


X  5[E  -E  -hAw  ]$  (Aw  )  +5 [E  -E+hAw  (-Aw  ) 

k'  k  cb  a  ab  k'  k  ad  b  ab 


P150960.PDF  [Page:  23  of  56] 

UNCLASSIFIED 

11 


+  <k  t  (g  +E  )  k'xk'  t  (c  +E  )  k> 
1  ac  a  k'  db  b  k' 


f  ©  (Am  )  $  (Am  )  ) 

cd  cd  cd  cd 


E  -E  -hw  E  -E  +hw 
k'  k  ac  k'  k  dby 


<k  fc  (e  +E  )  k'xk'  t  (g  +E  )  k> 
1  ac  c  k  1  db  d  k 


f  0  (Aw  )  9  (Aw  )  ) 

cd  ab  cd  ab 


E  -E  -hw  E  -E  +hw  . 

k '  k  ac  k 7  k  db^ 


+  <k  t  (c  +E  )  k'xk'  t  (g  +E  )  k> 

1  ac  a  k'  db  b  k' 


eciA"cd>  1 


E  -E  +hAw  E  -E  -hAw  . 
k'  k  ad  k '  k  cb' 


+  <k  t  (e  +E  )  k'xk'  t  (g  +E  )  k> 

ac  a  k  db  b  V 


f  0  (Aw  )  $  (Aw  ) 

cd  ab  cd  ab 


db  b  k 


E  -E  +hAw  E  -E  -hAw  .  , 

v  k '  k  ad  k '  k  cby  J 


[25] 


The  C  and  P  terms  of  eq.  25  (i*e.  the  last  four  terms)  still  contain 

off-shell  transition  matrices  (energy  argument  is  E  rather  than  E  )  as  well 

k'  k 

as  on-shell  matrices.  These  off-shell  transition  matrices  can  be  expressed  in 
terms  of  the  corresponding  on-shell  matrices.  The  two  are  connected  by  delay 
time  operators.  These  delay  times  are  well  known,  especially  through  the  work 
of  Felix  Smith  (Ref.  13).  They  are  present  in  all  collisions,  even  binary 
ones,  and  are  not  dependent  on  the  intervention  of  a  third  body.  The  delay 
times  express  the  fact  that,  in  a  binary  collision  for  example,  the  colliding 
pair  travels  for  a  time  as  a  quasi-bound  state  similar  to  a  dimer,  before 
finally  separating.  The  delay  time  is  a  measure  of  the  difference  in  time 
spent  in  collision  compared  with  the  time  over  the  same  trajectory  in  the 
absence  of  the  interaction.  It  can  be  positive  or  negative,  depending  on 
whether  the  interaction  is  attractive  or  repulsive.  It  is  therefore  to  be 
expected  that  this  type  of  kinematic  effect  will  appear  in  the  C  and  P  terms 
which  are  associated  with  collisional  overlap.  We  will  use  the  delay  time 
theory  to  account  for  this  effect  to  first  order  approximation. 


The  delay  time  At  is  defined  in  terms  of  the  energy  derivative  of  the 
phase  shift  7) 


At  =  2h  di)/d E.  [26] 

The  first  order  approximation  for  the  off-shell  t-matrix  products  of  eq.  25 
can  be  obtained  by  a  Taylor  series  expansion  about  the  on-shell  energies: 
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<k  t  (0  +E  )  k'xk'  t  (0  +E  )  k>  “  <k  t  (e  +E  )  k'xk'  t  (e  +E  )  k> 
1  aa  a  k'  1  1  bb  b  k'  1  1  aa  a  k  1  1  bb  b  k  1 


+  (E  -E  )  -5——  <k  t  (e  +E  )  k'>  <k'  fc  (e  +E  )  k> 
k'  k  SE  I  1  aa  a  k  1  I  1  bb  b  k  1 
k7  v  J 


+  (E  -E  )  <k  t  (0  +E  )  k'>  -3=* —  <k'  t  (c  +E  )  k> 
k'  k  1  aa  a  k  1  0Efc#  [  1  bb  b  k  1  J 


The  next  step  is  to  write  the  transition  matrix  elements  in  terms  of  the 
phase  shift.  The  phase  shift  is  complex  because  we  are  dealing  with  inelastic 
collisions  in  the  terms  C  and  P.  For  example,  assuming  spherical  symmetry  for 
simplicity,  and  dropping  the  internal  state  labels  for  the  same  reason, 
partial  wave  analysis  leads  to 

2i7) 

<k|t  (E  )  |k'>  =  -(2ir/m)  f(0)  (21+1)  (e  1-1)  P  (cos  0)  ,  [28] 

1  k'  1  mk  1 


<k|t  <Ek,)|k'>  =  -  ^  0(0)-  -  ^  0(0), 


where  the  wavenumber  k  =  k/h  and  we  have  defined  a  delay  time 


0(0)  =  2h  £  (21  +  1)  e  (07)  /SE,  )  P  (cos©)  . 

4-1  1  k'  1 


It  follows  from  eqs.  26-30  that,  in  general,  we  can  write  the  approximation 


<k  t  (0  +E  )  k'xk'  t  (0  +E  )  k> 

1  aa  a  k'  1  bb  b  k'  1 


^  <k  t  (e  +E  )  k'xk'  t  (0  +E  )  k> 

aa  a  k  1  1  bb  b  k  1 


2n  (E  -E  )  ,  *  * 

+  - - - —  Q(fi)/(Q)  +  £{Q)Qtn) 

2,  I  aa  bb  aa  bb 

m  k'  v 


It  is  now  convenient  to  express  all  the  t  matrix  elements  of  eq.  25  in  terms 
of  scattering  amplitudes  f(D,k')  and  introducing  eq.  31  we  obtain 


M(Aw)  =  n 
acdb  2 


2n  [ 
~  m  [ 


(Aw  )  f  (0)  -  *  (-Aw  J  fU 

a  ab  aa  b  ab  bb 
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2n 


m 


dk' 


4 


2ni  S[  E  -E  -hAw  1  $  (Aw  )+5[E  -E  +hAw  ]  3>  (-Aw  )  f(Q)  f  (Q) 

k7  k  cb  a  ab  k7  k  ad  b  ab  ac  db 


f  $  (Aw  )  +  ©  (Aw  )  #  (Aw  )  +  ©  (Aw  )") 

cd  ab  cd  cd  cd  cd  cd  ab 


E  -E  -hAw 

E  - 

v  k 7 

k  cb 

k' 

(E  f-E  ) 
k 7  k 

0  (Aw  ) 

cd  cd 

0  ^Aw 

2k' 

E  -E  -hAw 
k'  k  cb 

E  -E  ■ 
k 7  k 

ad 


f(Q)fJQ) 

ac  db 


$  (Aw  )  #  (Aw  )  "I 

cd  cd  cd  cd 


E  -E  +hAw  E  -E  +hW 
k 7  k  ad  k 7  k  dby 


x  Q(Q)f(Q)  +  jf  (ft)  Q  (ft) 
a  c  db  a  c  db 


AV 


[32] 


where  is  the  number  density  of  the  thermal  reservoir  molecules.  In 

passing,  we  note  the  difference  between  the  delay  time  Q  and  the  collision 
time  operator  T  (eq.  11)  .  Q,  as  defined,  relates  only  to  a  kinematic 
excursion  off  the  energy  shell  which  affects  only  the  translational  energy. 
It  is  not  directly  connected  with  the  intervention  of  a  photon  during  a 
binary  collision  process.  T,  on  the  other  hand,  is  defined  essentially  in 
terms  of  a  derivative  with  respect  to  the  whole  energy:  translational  plus 
internal  energy.  The  fact  that  the  off-shell  deviation  of  t(E+AE)  can  be 
expressed  in  terms  of  T  and  AE  =  hAw  (through  the  spectral  distribution 
functions  $)  is  due  to  the  direct  interaction  of  a  photon  with  the  absorber's 
internal  states  during  the  collision.  We  now  apply  the  optical  theorem  to  the 
forward  scattering  terms  of  eq.  32 


Im  f( 0)  =  k  cr l An  ,  [33] 

where  <r  is  the  total  cross-section.  Also,  remembering  that  k'  has  the 

same  magnitude  as  k,  but  is  at  angles  0  and  <p  to  it,  we  can  express  the 

integration  over  dk'in  eq.  32  as  one  over  dE  Thus, 

k' 

dk'  =  mk'  dE  d Q  /8n3,  [34] 

k' 


where  Q  is  a  solid  angle.  Equation  32  becomes 


M  (Aw)  =  n  i 
acdb  2 


2  IT  f 
-=-  [ 


$  (Aw  )Re  f  (0)  -  $  (-Aw  ) Re  f  (0) 

a  ab  aa  b  ab  bb 
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1 

+ — 
2m 


-*E 


$  (Aw  )<r  +  $ 

a  ab  aa  b 


dE  dQk' 
k ' 


i(i 


i  5[E  -E  -hAw 

*  k'  k  cb 


1 

2n 


(  V  (Aw  )  +  0  (Aw  ) 

ca  ab  cd  cd 


E  -E  -hAw 
k '  k  cb 


( E  -E  ) 
k '  k 

47Tk' 


r  eoia“cd> 


E  -E  -hAw 

k '  k  cb 


(-Aw 

ab 


]$  (Aw  )+5[E  -E  +hAw  ]$*(-Aw  )  f(Q)f*(Q) 

a  ab  k'  k  adb  ab  ac  db 


if  (Aw  )  +  0 

cd  cd  c 


E  -E  +hAw 
k'  k 


0  (Aw  )  if  (Aw  )  a  (Aw  )  " 

cd  cd  cd  cd  cd  cd 

-  -  -  +  - 

E  -E  -hw  E  -E  +hAw  E  -E  +hw 
k'  k  ac  k'  k  ad  k'  k  db"' 


f(n)fjn) 

ac  db 


x 


( 


Q(Q)fJQ)  +  f(Q)Q*AQ) 
a  c  db  a  c  db 


AV 


[35] 


Comparing  eq.  35  with  the  original  form  of  Fano's  operator  (eq.  2  with 
eqs.  3-8),  the  collision  time  theory  has  permitted  a  considerable  reduction 
in  the  number  of  terms.  The  terms  themselves  are  also  more  accessible  to 
physical  interpretation  and,  although  still  difficult  to  calculate  for 
specific  molecular  constituents,  their  effects  on  spectral  line  wings  can  be 
discussed -fairly  readily  in  analytical  terms. 


5 . 0  DIAGONAL  ELEMENTS  OF  THE  FANO  OPERATOR  AND  THE  LINE  WING 


For  an  analysis  of  the  non-Lorentzian  wings  of  isolated  spectral  lines 

we  are  primarily  concerned  with  the  diagonal  elements  M  (Aw)  .  The  analysis  of 

aabb 

the  non-diagonal  elements,  although  necessarily  more  complicated,  follows  the 
general  lines  of  the  following  treatment.  Denoting  Aw  by  Aw,  the  diagonal 

^  ab 

elements  are 


M  (Aw)  =  n 
aabb  2 


2n 

m 


$  (Aw) Re  f  (0)  -  $  (-Aw) Re  f(0) 

a  aa  b  bb 


--[■ 

2  L 

1 


$  (Aw)cr  +  $  (-Aw) 

a  aa  b 


- 1 

bbj 


2m 


dE  6Q  k' 
k' 
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2ft  ' 

hAw 

47rk' 


$  (Aw)  +  0  (Aw) 

;  ab  ab 


(  ©  (Aw) 

ab 


E  -E  -hAw 
•  k'  k 


E  -E  +hAw 
k'  k 


f(Q)f(Q) 

a  a  bb 


$  (Aw) 

ab 


E  -E  -hAw 
k '  k 


E  -E  +hAw 
k  '  k 


Q{Q)f(Q)  +  f{Q)Q(Q) 
a  a  bb  a  a  bb 


[36] 


J  y  AV 


This  result  shows  a  further  reduction  in  terms .  In  some  cases  it  may 

even  be  possible  to  neglect  the  last  term  (which  depends  basically  on  hAwQ) 

when  the  delay  times  are  very  small.  However,  the  most  interesting  aspect  of 

eq.  3  6  is  that,  because  of  the  factors  [E  -E  ±hAw]  ,  there  are  two  regimes 

k'  k 

for  the  diagonal  elements  M  (Aw)  ,  depending  on  whether  hAw  is  much  greater  or 

aabb 

2 

smaller  than  E  ~  k  /2m.  We  shall  first  examine  the  case  for  small  hAw  <<  E 
k  k 

which  is  the  near  wing  region. 


5 . 1  The  Near  Wing  Region 


The  near  wing  situation  was  discussed  in  Part  II  and  corresponds  to  the 
case  of  Aw  =  w  -  w  small  but  finite.  More  specifically,  when  hAw  « 

ab  ab  ab 

2 

k  /2m  “  kT  this  factor  can  be  neglected  simultaneously  in  both  the  delta 

functions  and  the  denominators  of  the  terms  in  eq.  36.  In  this  approximation 

the  principal  value  terms  cancel  out  completely  (as  indicated  by  Fano  in  Ref. 

1)  and  the  delay  time  term  is  negligible.  (In  fact,  the  justification  for  the 

dropping  of  all  these  terms  for  small  Aw  is  more  readily  evident  in  eq.  25)  . 

The  two  parts  of  the  singular  S  term  combine  under  a  single  delta  function 

5(E  -E  )  and  the  integration  over  E  can  be  carried  out.  Finally,  the 
k'  k  k' 

result  for  the  near  wing  version  of  eq.  36  can  be  written 


M  (Aw) 

aabb 


[ 


$  (Aw)  Ref  (0)  -  $  (-Aw)  Ref  (0) 

a  aa  b  bb 


XV 

~~2 


$  (Aw)  c 7  +  $  (-Aw) 

a  aa  b 


„  i 

bbj 


+  t*  (Aw  )  +$*  (-Aw  )  ]  an  f  to)  fjfi) 

2  a  ab  b  ab 


aa  bb  J 


1  . 


AV 


[37] 


where  v  is  the  average  perturber  velocity. 
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An  important  point  feature  of  this  result  is  that,  for  the  near  wing 

approximation,  the  frequency  dependence  of  M  (Aw)  is  determined  entirely  by 

aabb 

the  spectral  distribution  functions  $.  Although  we  require  a  knowledge  of  the 

pair  interaction  to  calculate  the  collision  time  T  (in  order  to  determine 

eq.  10)  ,  the  only  other  interaction  information  needed  is  at  the  impact 

theory  level  involving  f(0)  and  £{Q)  .  Thus  eq.  37  (with  eq.  10  for  $)  shows 

that  for  the  near  wing  region,  the  frequency  dependence  of  the  real  and 

imaginary  parts  of  M  (Aw)  are  determined  by  simple  rational  functions  whose 

aabb 

numerators  and  denominators  are  quadratic  in  Aw.  The  near  wing  theory  was 
applied  to  the  CO^  continuum  beyond  the  band  head  near  4  fim  in  Part  II. 

Because  the  range  of  validity  of  the  near  wing  approximation  is  about 
200  cm  (Aw  ^  kT/h)  for  ambient  temperatures,  xt  should  be  quite  accurate 
for  most  purposes  within  an  infrared  band  at  low  pressures.  In  such 
calculations  one  also  needs  to  take  into  account  the  line  coupling  effects  as 
described  in  I.  The  corresponding  frequency-dependent  non-diagonal  elements 
of  M  (Aw),  for  small  Aw,  were  given  in  Part  II. 

a  c  db 

5 . 2  The  Far  Wing  Region 


In  the  far  wing  region  where  hAw  »  E  we  cannot  make  the  same 

ab  k 

approximations  as  in  the  near  wing  region:  all  the  terms  of  eq.  36  must  be 

retained.  Although  the  spectral  distribution  functions  0  and  #  all  play  a 

role  in  the  frequency  dependence  of  M  (Aw)  ,  this  is  also  determined  by  the 

aabb 

specifics  of  the  perturber-absorber  interaction  which  affects  the  singular  as 
well  as  the  principal  value  terms.  However,  it  is  still  possible,  without 

going  into  the  details  of  the  interaction,  to  further  elucidate  the  nature  of 

these  terms  as  we  go  out  further  into  the  line  wing.  We  have  already 
indicated  that  the  nature  of  the  spectral  distribution  functions  0  and  $  is 
such  that  these  correspond  to  overlapping  binary  collisions  where  there  is  a 
mixing  of  the  absorber  internal  states  a  and  b. 

We  first  examine  the  effects  of  the  delta  functions  in  the  singular 

terms  of  eq.  36.  These  consist  of  S  and  S  where 
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dE  dfi  k'  8[E  -E  -hAw]  f  (£2)  f*  (Q)  ,  [38] 

k'  k'  k  a  a  bb 

dE  dQ  k'  S[E  -E  +hAw]  f{Q)  f  (Q)  .  [39] 

k'  k'  k  a  a  bb 

The  delta  functions  give  E  -»  (E  ±hAw)  .  Now,  in  the  partial  wave  analysis  of 

K  iC 

collision  processes  the  phase  shift  is  given  in  terms  of  the  energy  as  a 
function  of  the  wave  number  7c.  The  latter  is  related  to  the  linear  momentum  k 
by 

k  =  k/h  =  (2mE) 1/2 /h  . 

We  first  assume  that  Aw  is  positive  so  that  we  are  looking  at  the  high 

2 

frequency  wing.  If  hAw  »  k  /2m  then  E  -»  (E  ±hAw)  means 

k'  k 

2  1/2 

*'  •*  [  2m  (k  /2m  +  hAw  )  ]  /h 

and 

2  1/2 

7c '  [  2m  (k  /2m  -  hAw  )  ]  7  /h 

Thus  in  the  wave  number  increases  with  Aw  and  the  photon  energy  hAw  is 
absorbed  into  the  translational  energy.  In  S  the  wave  number  7c'  becomes 
imaginary . 

Certain  conclusions  may  now  be  drawn  on  the  nature  and  relative 

contributions  of  the  terms  S+  and  S  in  the  far  wing  region.  In  the  near  wing 

region  these  terms  merge  together  into  a  single  term  where  the  dependence  on 

* 

Aw  is  given  solely  through  the  function  [$  (Aw)+$  (-Aw)].  On  passing  into  the 

a  b 

far  wing  region  the  nature  of  the  two  terms  becomes  quite  different.  Term  S 

+ 

continues  to  behave  as  a  purely  scattering  function  with  scattering 
amplitudes  corresponding  to  continuously  increasing  wave  number:  the  photon 
energy  hAw  passing  into  the  translational  energy.  In  the  term  S  we  have  to 
remember  that  we  are  dealing  with  elastic  collisions  and  therefore  the  phase 
shifts  are  basically  real.  (The  situation  is  more  complicated  for  the 
corresponding  part  of  the  non-diagonal  elements  M.(Aw)  which  concerns 

acdb 

inelastic  collisions  and  the  phase  shifts  are  complex) .  With  increasing  Aw, 


h~ 1  (2mhAw) 1/2 


+  2mhAw] 


1/2 


[40a] 


.  -1  ^  »  1/2 
ih  (2mhAw) 


(i _ hl_) 

[  2mhA(t>J 


1/2 


[40b] 


=  —  $  (-Aw) 

2m 


S  =  —  $  (Aw) 

2m  a 
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in  the  transition  from  the  near  wing  to  the  far  wing,  the  delta  function 

decreases  the  wave  number  and  puts  the  collision  into  a  low  energy  regime 

21+1 

where  the  phase  shifts  17  are  proportional  (Ref.  14)  to  k  .  When  hAw  = 
2 

k  /2m,  k'  =  0  and  the  phase  shift  goes  to  zero.  At  this  point,  S  is  zero. 

Beyond  this  point  the  wave  number  and  phase  shift  are  imaginary.  The 

appearance  of  imaginary  phase  shifts  signifies  a  transition  from  a  scattering 

process  to  one  of  absorption.  For  these  reasons  we  can  conclude  that  the 

functional  dependence  of  M  (Aw)  on  Aw  changes  dramatically  in  the  transition 

aabb 

region  which  will  be  marked  by  a  turning  point  in  the  wing  shape.  A  further 

conclusion  is  that  the  far  line  wings  are  normally  asymmetric.  This  is 

because  on  the  low  frequency  side,  where  Aw  is  now  negative,  the  roles  of  S+ 

and  S  are  interchanged:  S  corresponding  to  absorption  and  S  to  scattering. 

-  +  - 

The  analysis  of  the  overlapping  collision  terms  is  more  difficult  in  the 

absence  of  specific  calculations.  We  can  see  however,  that  the  two  poles 

given  by  [E  -E  ±hAw]  1  which  were  near  coincident  in  the  near  wing  and  thus 
k'  k 

cancelled  out,  now  move  in  opposite  directions:  one  towards  infinite  energy 
and  the  other  towards  zero  energy.  There  will  therefore  be  an  increasing 
contribution  from  these  terms  as  we  move  into  the  far  wing  region.  The 
relative  importance  of  this  contribution  is  not  evident. 


Having  established  that  there  is  a  turning  point  near  the  transition 
region  between  the  near  and  far  wings  of  an  isolated  line,  it  would  be  useful 
to  determine  the  direction  of  this  feature,  i.e.,  whether  the  far  wing  dips 
downwards  or  simply  flattens  out.  An  indication  of  this  can  found  from  the 
asymptotic  limit  of  M  (Aw)  as  Aw  00.  This  limit  can  be  determined  from  eq. 

aabb 

36  as  follows.  First  of  all  we  take  the  asymptotic  limits  of  $(Aw) ,  ©(Aw)  and 
${Aw),  which  are  constants  from  eqs .  10c,  lOd  and  21,  and  neglect  E^  compared 
with  hAw.  Equation  36  then  gives: 


Lt  M(Aw) 

.  aabb 

Aw**  co 


-  V 


271 


$  (CO) Kef  (0)-$  (-co)Kef  (0) 

a  aa  b  bb 


]-*[■ 


$  (co)cr  +$  (-00) 
a  aa  b 


,  1 

bbj 


2m 


dEdQ  k  5[E-hAw]$  (co)  +8  [E+hAw]  $  ( 
a  b 


— co)  j 


£{Q)f(Q) 

a  a  bb 
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hAw 


0(«) 


Tim 


dEdfl  k 


_2  T-2*  2 

E  -h  Aw 


f  (Q)  fJD) 

a  a  bb 


hAw 

-  © (oo) 

27rm 


- 

P  1 

, 

dE  E  dQ 

T-,2  ,2a  2 

Q(Q)f*(£2)  +  f(£2)Q*(£2) 

a  a  bb  a  a  bb  J 

J 

E  -h  Aw 

^  J 

-  . [41] 

AV 


where  the  integrations  over  E  are  from  0  to  oo.  The  first  two  terms  are 
constants.  To  evaluate  the  asymptotic  limits  of  the  remaining  terms  we  can 
use  approximations  for  the  scattering  amplitudes  f(Q,k)  at  low  energies  and 
high  energies.  The  Born  approximation  (Ref.  14)  can  be  used  for  the  phase 
shift  for  large  E  {=  k  /2m  =  Jc2/2mh2)  -*  co.  The  first  Born  approximation  gives 
f(Qrk)  ~  k  and  the  third  term  has  components  which  fall  off  as  Aw  4  and 
therefore  gives  zero  at  Aw  =  co.  Higher  order  Born  approximations  fall  off 
even  faster.  The  fourth  and  fifth  terms  can  be  written  as 


I 

4 


hAw  r 

Lt  -  0  (oo) 

Aw-*oo  7rm 

J  0 


dE  E 


1/2 


„2  ,2fl  2 

E  -h  Aw 


F(\/E) 


I 

5 


hAw 

Lt  -  ©  (oo) 

Aw->oo  2  Tim 


03 

dE  E 

J  0 


[E2-h2AW2J 


G(/E) 


At  the  upper  limit,  the  functions  F^E)  and  G(/E)  can  be  approximated  using 

the  first  Born  approximation  where  both  converge.  At  the  lower  limit  we  can 

21+1 

use  a  low  energy  approximation  (Ref.  14)  for  the  phase  shift  7}  ~  k  .  With 

/  41  ,  41-1 

7?1  small,  as  k  ->  0,  we  have  terms  in  1  with  F(vE)  ~  k  and  G(vE)  -  k 
Again,  both  integrals  are  convergent  at  the  lower  limit. 


There  remains  the  problem  of  the  pole  at  E  =  hAw  in  each  integral.  Since 

both  integrals  fall  off  rapidly  in  the  regions  of  low  and  high  energy  we  can 

expect  that  the  main  contribution  to  the  integrals  I  and  I  will  in  fact 

4  5 

come  from  the  pole.  However,  since  we  are  only  interested  in  the  limiting 
case  hAw  very  large  it  is  possible  to  approximate  the  functions  F(\/e)  and 
G(,/e)  by  their  values  at  E  =  hAw  and  take  them  out  of  the  integration.  More 
exactly,  we  develop  the  functions  about  this  value  using  a  Taylor  series 
expansion 


F(v'E)  =  F(v'hAu)  +  (vWhA^)F'  (VE5S)  +  (v'E-VhAw)2F"  (VhAw) /2  +  .  .  .. 
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Then,  provided  that  the  first  and  second  derivatives  exist  and  fall  off  with 

hAw  faster  than  the  F(V'hAw)  we  get  an  accurate  approximation  to  the  integral. 

These  conditions  are  valid  for  both  functions  F  (-/hAw)  and  G  (VhAw)  using  the 

_  _2 

first  Born  and  higher  approximations  since,  for  example,  F(VhAw)  ~  (hAw) 

The  principal  values  at  the  poles  can  now  be  evaluated  as 


hAw 


~E  +  e 


I  -  -  Lt  -  ©  (co)  F  ("/hAw)  P  ,  with  P  =  Lt 

4  A  4  4 

Aw+oo  7rm  E**oo 


dE  E 


1/2 


_  2  ,  2  .  2 
E  -h  Aw 


hAw 


~E  +  e 


I5  “  -  Lt  -  0  (oo)  G(\/hAw)  P  ,  with  P  =  Lt 

Aw-*oo  2  Tim  e-»oo 


dE  E 


E-e 


_2  ,  2  A  2 
E  -h  Aw 


and  the  limit  e  ^  0  is  to  be  taken  after  the  limit  E  *»  oo.  We  obtain 


P  =  Lt 
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E-*co  V'hAw 
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-  log 
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vhAw, 
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E+oo  2 


-  |  log | (E+e) 2  -  (hAw) 2 |  -  log|(E-e)2  -  (hAw)2|  j  = 


The  fourth  and  fifth  terms  are  therefore  also  zero  and  eq.  (41)  becomes 

y  y* 
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r  j  u 

-  —  — -  Re  f  (0)  -  —  Re  f 

m  L  5,  aa  y*  ] 


(0) 

bb 


IV 

2 


r  y  y  *i  \ 

+  -1  <r  ]\  . 

U  aa  T*  bbJiAV 

a  b 


[42] 


where  we  have  used  the  asymptotic  expressions  of  eqs .  10c  and  lOd.  Thus,  the 

real  and  imaginary  parts  of  M  (co)  converge  to  constants  which  are  simple 

aabb 

functions  of  the  elementary  collisional  constituents  f(0)  and  the  collision 
times.  A  good  approximation  to  the  right-hand  side  of  eg.  42  can  be  easily 
obtained.  Thus,  assuming  that  all  the  parameters  in  the  two  states  are 
approximately  equal,  the  first  term  is  small  and  can  be  neglected  and  we  have 
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assuming  that  c r  “  cr  ^  <r.  The  average  of  n  vcr  can  be  identified  with  the 
aa  bb  2 

collision  frequency  T^f  of  the  kinetic  theory  of  gases,  giving 


* 


a  b 


[42a] 


The  theory  predicts  that  after  the  turning  point  near  the  transition 

region  the  real  and  parts  of  M  (oo)  tend  towards  constant  values.  This  results 

aabb 

in  a  Voigt  profile  for  the  far  line  wing.  To  see  this  we  write 


M  (w)  =  A  (w)  -ir  (0)) 

aabb  aabb  aabb 


and  the  absorption  coefficient  for  an  isolated  infrared  line  (of  Doppler 
width  oM  in  the  far  wing  is  given  by 

- 1  r  oo 

-3/2 

a(oj)  ~  n  n  S  w  tanh(fthw/2)  dt  e  -  ,  [43] 

1  [A  0)  -  A  (oo)  -  o'  t]2+  r  (co)  2 

-oo  D 


r  y 

y  i 

r  y 

yfa  1 
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T  J 
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T  J 
b 

Some  experimental  corroboration  is  available  for  this  result.  In  Part  II 

we  applied  the  near  wing  theory  to  the  continuum  beyond  the  band  head  of  the 

v  band  of  CO^  near  4.2  fxm.  As  discussed  in  Part  II  the  experimental  results 

of  Cann  et  al.  (Ref.  15)  were  used  for  CO^  broadened  by  at  a  temperature  of 

29  6K.  Reference  15  also  gives  reults  for  O  and  self -broadening .  These 

^  -1 

results  are  in  the  form  of  a  form  factor  ^(Ao>)  which,  for  Ao>  >  50  cm  ,  is 
basically  equivalent  to  T  (to)  (as  discussed  in  Part  II).  Cann  et  al . 
determined  £(Aw)  empirically  by  an  inversion  of  the  continuum  absorption 
assuming  that  the  same  ^(Aw)  was  valid  for  each  line  in  the  CO^  band.  This 
form  factor  is  plotted  in  Fig.  1  for  self  and  broadening  (O^  and  N2  are 
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FIGURE  1  -  Form  factors  for  self-  and  foreign-gas  broadening  of  CO^  (Ref.  15) 

near  identical)  .  In  each  case,  the  predicted  turning  point  is  seen  near  ~ 
200  cm  1  and  the  ^(Aw)  then  levels  off  towards  a  constant  value.  In  the  case 
of  CO  broadened  by  N  ,  this  value  is  about  0.004.  This  result  can  be  used  in 
conjunction  with  those  of  Part  II  where,  applying  the  near  wing  approximation 
eq.  37  to  the  experimental  data  of  Cann  et  al .  ,  estimates  were  obtained  for 
the  real  and  imaginary  parts  of  T  and  T  by  inversion  (assuming  that  *3  << 

aa  bb 

T)  .  We  obtained  Re  ( T  ,T  )  “  24,  27  x  10  s  and  Im  ( T  ,T  )  ^8,  7  x 

aa  bb  aa  bb 

10  13s.  Taking  -  0.04  cm1  and  an  average  half  width  r^O^  “  0-07  cm  , 
and  assuming  the  same  relative  proportions  for  the  real  and  imaginary  parts 
of  S'  as  for  T,  we  find  that  the  free  propagation  times  are  smaller  than  the 
collision  times  T  by  a  factor  ^  7  x  10  3 .  This  result  also  justifies  the 
approximation  of  Part  II  that,  in  the  near  wing  approximation,  we  can  neglect 
S'  in  comparison  with  T  in  the  spectral  distribution  functions  $. 


6.0  DISCUSSION  AND  CONCLUSIONS 


The  collision  time  theory  allows  a  considerable  simplification  of  the 
number  and  form  of  the  terms  contributing  to  the  matrix  elements  of  Fano's 
relaxation  operator,  particularly  for  the  diagonal  elements  T^€ 
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analysis  of  the  diagonal  terms  and  their  application  to  an  isolated  line  wing 
shows  that  the  latter  is  composed  of  two  regions  of  different  functional 
behaviour:  the  near  and  far  wing  regions.  The  demarcation  between  the  two 
regions  occurs  near  a  displacement  from  the  line  centre  Aw  ~  kT/h.  In  the 
near  wing  the  collision  parameters  are  essentially  those  of  the  impact 
approximation  and  the  shape  of  this  part  of  the  wing  is  determined  by  the 
rapidly  changing  spectral  distribution  functions  $  (Aw)  and  $  (-Aw)  .  Near  the 

a  b 

transition  region/  between  the  near  and  far  wings,  several  changes  occur.  The 

singular  contributions,  S+  and  S  ,  which  were  merged  together  in  the  near 

wing,  begin  to  behave  divergently:  continuing  to  .behave  as  a  scattering 

function  and  becoming  an  absorption  function.  At  the  same  time,  the 

collision  overlap  terms  cease  to  cancel  out  mutually  and  become 

non-negligible .  The  shape  of  the  far  wing  is  largely  determined  by  the  nature 

and  specific  details  of  the  binary  molecular  interactions.  This  is  because 

■* 

although  the  spectral  distribution  functions  $  (Aw),  $  (-Aw),  0(Aw)  and  #{Aw) 

a  b 

exhibit  significant  variations  at  small  values  of  Aw,  in  the  far  wing  they 
tend  to  level  out  to  constant  values. 

All  theories  of  line  shape  must  obey  the  Nyquist  relation 

‘°°  X"  («> 

d0>  -  =  a  £  S  , 

tanh(/3hw/2)  ab 

CO 

where  X"  (w)  is  the  imaginary  part  of  the  susceptibility  and  S  is  the  line 

ab 

strength.  This  condition  is  very  important  because  it  ensures  that  the 

integrated  band  intensity  is  unaltered  by  a  redistribution  of  line  intensity, 

due  to  line  coupling  effects  and  strong  deviations  of  individual  line  wings 

from  the  Lorentz  form.  Now,  it  does  not  appear  possible  to  prove  with 

*  complete  mathematical  rigour  that  the  theory  obeys  the  Nyquist  relation.  This 

is  because  the  expressions  for  the  frequency  dependence  of  the  relaxation 

matrix  M(Aw)  are  too  complicated  for  analytical  evaluation  in  the 
acdb 

integration.  However,  it  is  fairly  easy,  using  asymptotic  properties  that  we 

have  derived  for  M(Aw)  ,  to  carry  out  a  good  analytic  approximation  of  the 

acdb 

Nyquist  integral  to  show  that  it  tends  towards  the  desired  result.  This 
analysis  is  carried  out  in  Appendix  B,  where  it  is  evident  that  the  main 
reason  that  integral  tends  towards  the  Nyquist  condition  is  due  to  the  nature 


1 

2 
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of  the  line  coupling  mechanism  derived  in  Part  I. 

Another  test  of  the  validity  of  the  theory  is  the  proper  convergence  of 
* 

moment  integrals  .  Thus,  the  moments  of  the  real  part  of  the  complex 
susceptibility  must  be  finite  and  the  question  arises  whether  this  is 

true  in  the  present  case,  where  the  asymptotic  limits  of  the  line  wing  behave 
as  Voigt  profiles,  eq.  43.  This  is  demonstrated  in  Appendix  C. 

Whilst  the  collision  time  theory  provides  a  considerable  reduction  of 
Fano's  relaxation  matrix  and  enables  a  better  understanding  of  the 
significance  of  the  various  terms,  entirely  theoretical  calculations  are 
still  extremely  difficult,  especially  for  molecular  interactions.  For  the 
near  wing,  assuming  impact  parameters  are  already  available,  it  is  necessary 
to  calculate  the  collision  times  T  and  S'  (diagonal  elements) .  In  the  far  wing 
we  also  have  to  perform  new  types  of  calculations  to  take  into  account  the 
collisional  absorption  and  overlapping  binary  effects. 

In  terms  of  extracting  collision  interaction  information  from 
experimental  line  wing  results  by  inversion,  the  near  wing  yields  estimates 
only  of  the  collision  time  T.  Information  about  the  relative  magnitude  of  the 
free  propagation  time  S',  compared  with  T,  can  be  obtained  from  the  asymptotic 
far  line  wing. 

In  the  absence  of  specific  calculations,  direct  comparison  with  other 
theories  is  not  yet  possible.  A  comparison  with  Rosenkranz's  quasistatic 
theory  (Ref.  8)  and  its  further  development  (Refs.  9  and  10)  would  be 
particularly  interesting.  This  is  because  it  applies  to  the  same  far  wing 
region  but  is  based  on  diametrically  opposed  concepts  and  assumptions.  Thus, 
a  basic  premise  of  quasistatic  theories  is  that  the  collision  times  are 
infinite  and  the  statistical  distribution  can  be  treated  as  though  the 
molecules  were  stationary.  The  theory  is  therefore 

primarily  applicable  to  the  asymptotic  far  wing.  However,  very  recently  Ma 
and  Tipping  (Ref.  16)  have  further  developed  the  quasistatic  theory  to 

include  corrections  for  the  near  wing  and  molecular  motion. 

* 
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In  contrast,  the  collision  time  theory  used  in  this  report  and  in  Part 

II  is  very  different  in  concept  to  the  quasistatic  theory.  The  key  issue  is 

that  the  theory  is  developed  strictly  in  the  frequency  domain  where  the  Fano 

relaxation  operator  is  defined.  In  this  domain,  "times"  appear  as  constants 

or  parameters  and  not  variables.  The  principal  variable  is  frequency  in  the 

form  of  a  quantum  of  energy  hAo).  This  is  borne  out  by  the  collision  time 

theory  which  relates  deviations  off  the  energy  shell  of  the  transition  matrix 

to  collision  times  which  are  fixed  calculable  quantities.  As  Fano  noted  (Ref. 

1)  the  variation  of  the  transition  operator  t(E)  as  a  function  of  its  energy 

argument  E  is  a  measure  of  the  duration  of  a  collision.  The  collision  times 

which  arise  out  of  the  collision  time  theory  do  so  through  a  quantum 

mechanical  time  operator  -ihS/dE  (see  II)  .  The  collision  duration  time  T  is 

therefore  a  quantum  mechanical  quantity  with  the  dimensions  of  time.  Its 

matrix  elements  do  not  have  an  exact  classical  counterpart:  T  is,  in  fact, 

complex.  Similarly,  when  we  refer  to  the  free  propagator  time  operator  3* 

(defined  by  eq.  12)  as  the  time  a  molecular  pair  would  spend  in  the  same 

vicinity  in  the  absence  of  the  interaction,  this  is  only  true  in  a  quantum 

mechanical  sense.  S'  is,  in  fact,  merely  a  convenient  form  in  which  to  express 

the  free  propagator  (E-H  )  1:  it  has  no  classical  counterpart  and,  like  T,  is 

o 

complex . 
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APPENDIX  A 


Evaluation  of  some  Principal  Value  Integrals 


Consider  first  the  integral  I  : 
„oo 

P 


nil  = 
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dip  - 

ip-e  -E  , 

oo  c  k ' 
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(  i  (T  /h)  (ip-e  -E 
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v  H  c  k  '  J 


I  l-i  (T  /h)  (ip-e  -E  )  I 
v  c  c  k'  J 

Defining  the  complex  parameters 

A  =  T  /T  ,  B  =  h/T  ,  C  =  T  /T* ,  D  =  h/T* , 
c  c  c  d  d  d 

and  putting  A  =  hAo)  ,  x  =  c  +E  ,  I  can  be  written 
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where  we  have  used 


dip 


=  0 .  Changing  the  variables  in  the  remaining 


i/f-x 


integrals  (which  are  not  principal  values) 
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Now,  B  -  a+ib  and  D  -  c-id,  whence  integrals  are,  for  Re  p  >  0,  of  form 
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Whence,  eg.  A3  becomes 


I  =  A  ( 1-C )  +- 


iACD  iCD  iACD  1  2iACD  iCD 

-  + - 4 -  =  A(l-C)  + - .  [A4] 
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However,  the  last  term  of  eq.  A4  is  related  to  $  (-Aw  )  , 
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X  =  (1+A)  (1-C)  + - $*(-Aw  ). 
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It  follows  that  I  is  the  same  as  I  with  hAw  replaced  by  hAw  . 

2  1  cd  ab 

In  an  exactly  similar  fashion  I  and  I  follow.  Thus,  we  obtain: 
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It  is.  now  convenient  to  define  two  new  spectral  distributions: 
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APPENDIX  B 


Approximate  Evaluation  of  the  Nyquist  Relation 
We  wish  to  show  that  the  theory  conforms  to  the  Nyquist  relation 


X"  («) 

dw  -  =  n  £  S  . 

tanh(/3hw/2)  1  ab  ab 


In  Part  I  it  was  shown  that  the  integral  on  the  LHS  is  basically 
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where,  for  simplicity,  we  have  neglected  the  Doppler  width.  This 
approximation  does  not  effect  the  discussion  and  can  be  readily  omitted.  The 
function  F(a>)  is  given  by  eq.  38  of  Part  I: 
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(Note  that  in  eqs .  B3  -  B4  we  have  retained  the  notation  of  Part  I  for  the 
arguments  of  A  and  T,  i.e.  A(o>)  rather  than  A(Afc>)  which  was  adopted  in  the 
present  paper).  P(w)  and  Q(a>),  and  the  H(w)  and  J(o>)  ,  are  coupling 

ab  ab  acdb  acdb 

"coefficients"  and  would  be  constants  were  it  not  for  their  frequency 

dependence  on  such  quantities  as  A(±w)  and  ]H(±a))  etc.  We  consider  the  first 

aabb  aabb 

term  in  eq.  B2 ,  i.e.,  the  positive  resonance  term  F(w)  .  The  negative 

resonance  term  can  be  evaluated  in  an  identical  manner.  Making  use  of  the 
relations  between  the  coupling  coefficients,  namely 


P(p> 

ab 


E  '  and 


c^a 

d=£b 


Q(w) 

ab 


=  E 


c^a 

d^b 


J(0>1 

acc 


it  is  convenient  to  write  the  integral  in  the  form 
1 

I  =  n  -  y  S 

1  1  n  ^  ab 

2  ab 

-oo 

da)  -( 


do)  F(u>) 
ab 


=  n  —  £  S 

2n  ab  a 


ri&> 


I  [w  -  w  -A(w)  ]2  +  r2(«) 

v  ab  aabb  aabb 


+  E 


c^a  . 

d*b  L 


ri&, 


ri!&. 


I  [w  -  u  -A(o>)  ]2  +  r2(w)  [a  -  0)  -A(w)  3 2  +  r2(w)J 

v  ab  aabb  aabb  cd  ccad  ccda/ 


+  J(w) 

acdb 


[w-w  -A  ] 

ab  aabb 


[w-w  -A(w)  ] 

cd  ccdd 


I  [W  "  Wab-A i^b]  +  L Ifi,-."cd-A^d3  +  J 


[B8] 


First  of  all,  from  the  asymptotic  properties  of  the  diagonal  and  non-diagonal 
elements  of  A(a>)  and  T(w)  derived  in  this  paper  as  0)  oo,  it  is  evident  that 
the  integral  converges.  Thus  the  diagonal  elements  tend  to  small  constants 
and  the  non-diagonal  elements  tend  to  zero . 


In  order  to  carry  out  an  approximate  evaluation  of  the  integral  the  only 
latitude  we  have  is  to  replace  the  diagonal  elements  and  and  the 

non-diagonal  elements  T  (w)  and  A(w)  by  constant  values  such  that  I  is 

acdb  acdb  i 

dominated  by  the  contributions  at  these  values.  In  fact,  as  we  shall  see,  we 

need  only  to  make  appropriate  choices  for  the  diagonal  elements  and 

A  (0))  .  Consider  first  the  first  term  of  eq.  B8 .  From  the  present  paper  and 
aabb 
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Part  II  we  know  that  r(w)  has  a  peak  value  at  the  line  centre  which  is  its 

aabb 

impact  approximation  limit  HO)  ,  so  this  is  the  value  we  choose.  It  can 

aabb 

therefore  be  taken  out  of  the  integral.  A(w)  ,  on  the  other  hand,  is  small  at 

aabb 

the  line  centre,  but  exhibits  a  maximum  and  a  minimum  on  either  side  (see 

Part  II)  .  The  contributions  from  the  maximum  and  minimum  balance  out  under 

the  integration,  so  it  is  a  good  approximation  to  replace  A ( Ct>)  by  A(0)  .  The 

aabb  aabb 

same  approximations  are  also  made  for  the  same  functional  components  in  the 
remaining  terms  of  eq.  B8 .  The  A(u>)  and  T  {0))  may  manifest  satellite  line 
maxima  arising  from  any  low  energy  resonances  in  the  scattering  amplitudes, 
however  these  secondary  maxima,  if  they  occur,  must  be  very  small  and  we  can 
assume  that  they  contribute  little  to  the  integrand. 


There  remains  the  question  of  determining  the  appropriate  approximate 

values  for  the  coefficients  H(&>)  and  J(w)  .  Since  these  depend  entirely  on 

acdb  acab 

all  the  diagonal  and  non-diagonal  elements  through  A(tt)  and  B(o>)  of  eqs .  B6 

acab  acab 

and  B7  a  reasonable  choice  would  be  to  use  the  impact  approximation  values  to 

obtain  values  H(0)  and  J(0)  .  Substituting  these  values  in  eq.  B8  we  obtain 
acdb  acdb 

an  integral  which  can  be  readily  evaluated: 
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+  jmb 


[W-W  -A (0}  ] 
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l  [03-03  -A ( o )  ] 2  +  r2(0)  [«  -  w  -A(0)  ] 2  +  r2(0)  l 

v  ab  aabb  aabb  cd  ccad  ccda  7 
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2ti  ab  V  c^a  L  V  }  ^  2a 

d^b 

The  same  result  is  obtained  from  the  negative  resonance  term  and  the 

Nyquist  is  therefore  valid  to  a  good  approximation.  Equation  B9  shows  that  it 

matters  little  what  "reasonable"  approximation  is  used  for  H(w)  and  J  (o))  , 

acab  acab 

since  the  terms  associated  with  these  cancel  out.  This  is  due  to  the  form  of 


the  line  coupling  mechanism  which  causes  the  gain  and  loss  coupling  terms 
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tend  to  cancel  out  over  the  overall  band  intensity. 


APPENDIX  C 


Moments  of  the  Susceptibility  in  the  Far  Wing  Asymptotic  Limit 


We  wish  to  show  that  the  moments  of  the  real  part  of  the  complex 
susceptibility  =  xf  (k>)  -i x"  (w)  ,  are  finite.  The  response  function  <f>  (t) 

and  X(u)  are  related  by  the  following  expression  (see  Part  I,  eq.  20) 

r00 

x(u>)  =  I  dT  e"lWt  <p  (t),  [Cl] 

O 

from  which,  by  inversion  and  successive  time  derivatives  (see  for  example 
Kubo,  Ref. 17),  the  moments  are  defined  by 


(-l)m  0(2m)  (0)  = 

W1 


_  2m  „  ,  . 

do)  o)  x  (a>)  . 


The  imaginary  part  of  the  susceptibility  (see  eqs .  34  -  38  of  Part  II)  is 


n  c  tanhOho)/2)  f00  -t2 

*"  l">  -  - — —  e  s.b  dt  -  +  1  • 

Sn  ab 

-oo 


Retaining  only  the  positive  resonance  term  F(o>)  ,  and  in  particular  the 
isolated  line  (neglecting  line  coupling  effects) ,  we  are  concerned,  in  the 
evaluation  of  the  moments  Eq.  C2 ,  with  integrals  of  the  form: 


r00  r“  -t2  r<w) 

J  =  -  dw  w  tanh  (|3hw/2 )  dt  e  - - - - - 

m  71  [w-w  -A(W)-<r  t]2  +  T  (0)} 

-00  -  00  o 


dw  w  tanhOhw/2)  U(x,y). 


where  U(x,y)  is  the  real  part  of  the  complex  error  function  and  the  variables 


=  [«  -  w  -  A(w)]/cr  ,  y  =  T (w)  /c  ,  v  =  Doppler  width. 

o  D  D  D 
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The  problem  is  to  show  that  the  integrand  (which  is  an  odd  function  due  to 
the  tanhOhw/2)  function  which  is  odd)  of  eq.  C4  is  finite  at  the  limits  10  -» 
±oo  where,  as  we  have  shown,  r(co)  and  A(a>)  tend  to  constants  T  and  A  in  the 
asymptotic  far  wing.  This  is  readily  evident  using  an  asymptotic  expansion 
for  the  function  U(x,y)  derived  by  Reichel,  Ref.  18: 


U  (x,y) 


n  =  o  n  I 


S  , 
n 


The  coefficients  S^(y)  axe  functions  of  y  only  and  obey  the  recurrence 
relation 


S 

n 


(2n-l)  (1  Sn-i)' 


where  S 

o 


1/2 

n  y 


erfc  y. 


In  the  limit  U),  x  ■*  ±oo  the  integrand  of  eq.  C4 


2m  2  2n 

2m  X  -x  *  x 

U)  tanhOhw/2)  U(x,y)  ->  ±  -  e  T  -  S  ->  0, 

1/2  ^  ,  n 

TC  y  n  =  o  n! 

because  of  the  exponential  factor.  The  line  coupling  terms  in  F(o>)  contain 
the  imaginary  part  V(x,y)  of  the  complex  error  function,  and  one  can  use 
Reichel ' s  form 


V  (x,  y) 


x 


1/2  2 
n  y 


2 

-x  oo 


2n 
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z  — 

n  =  o  n! 
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n  +  l 


to  show  that  their  contribution  to  the  moments  is  also  finite. 
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